- 1 - 



FJ-F263 



APPARATUS FOR CALCULATING IMMUNITY FROM RADIATED 
ELECTROMAGNETIC FIELD. METHOD FOR ACHIEVING CALCULATION, 
AND STORAGE MEDIUM STORING PROGRAMS THEREFOR 



BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an apparatus 
and method for calculating immunity from a radiated 
electromagnetic field which make possible high speed 
simulation of an electric current flowing through an 
electronic apparatus due to a radio wave radiated from an 
antenna and to a storage medium storing programs used for 
the same. 

One of the new requirements being imposed by 
society on electronic apparatuses is that they not be 
affected by radio waves of under a certain level radiated 
from other electronic apparatuses . Tough regulations on 
this have been established in the major countries of the 
world . 

Under radio wave regulations , examinations are 
conducted to determine if electronic apparatuses are 
affected by radio waves radiated from antennas. This has 
made necessary the development of technology for 
simulation of the effect of radio waves radiated from 
antennas on electronic apparatuses . 

2 . Description of the Related Art 

The electric current and magnetic current 
flowing through parts of an object can be theoretically 
found by solving Maxwell's electromagnetic equations 
under given boundary conditions. 

As a method of solving this, there is the 
moment method. The moment method is one of the methods of 
solving integration equations derived from Maxwell's 
electromagnetic equations by segmenting an object into 
small elements and therefore is able to handle any three- 
dimensionally shaped object. As a reference on the moment 



method, there is "H. N. Wang, J. H. Richmond, and M. C. 
Gilreath: "Sinusoidal reaction formulation for radiation 
and scattering from conducting surface", IEEE 
TRANSACTIONS ANTENNAS PROPAGATION, vol. AP-23, 1975". 

The electric current and magnetic current 
flowing through elements are found by segmenting the 
configuration of an apparatus to be simulated into 
meshes, selecting a frequency to be processed, finding 
the mutual impedance, mutual admittance, and mutual 
reaction among the mesh-like elements for the selected 
frequency by performing predetermined computations, 
substituting the found mutual impedance etc. and a wave 
source specified by the configuration information into 
simultaneous equations under the moment method, and 
solving those equations . 

That is, when handling a metal object, the 
metal portion is segmented into meshes as the object of 
analysis, the mutual impedance Z (value at frequency 
being processed) among mesh-like metal elements is found, 
and the simultaneous equations of the moment method 
standing among the mutual impedance Z Lj , a wave source V ± 
of that frequency component, and an electric current I t 
of that frequency component flowing through the meshed 
metal elements: 

[ZijlIIil = [Vi] 

where, [] indicates a matrix 
are solved to find the electric current I L flowing 
through the metal elements . 

Note that the mutual impedance shows the 
relationship between the electric field induced by the 
electric current flowing through an element and the 
electric current flowing through another element. The 
mutual admittance becomes necessary when considering the 
existence of a dielectric and shows the relationship 
between a magnetic field induced by a magnetic current 
passing through one element and the magnetic current 
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passing through another element. The mutual reaction 
becomes necessary when considering the existence of a 
dielectric and shows the relationship between the 
electric field (magnetic field) induced by an electric 
5 current (magnetic current) flowing through one element 

and the magnetic current (electric current) passing 
through another element. An electric current flows 
through metal, while an electric current and magnetic 
current flow on the surface of the dielectric . 
10 Up until the present time, the fact is that no 

technology for simulation of the effect of a radio wave 
radiated from an antenna on an electronic apparatus had 
:p been developed. 

It has however now become possible to simulate 
15 the effect of a radio wave radiated from an antenna on an 

electronic apparatus by using the moment method. 

That is, the electronic apparatus to be 
simulated and the antenna radiating the radio wave are 
set as a single object for application of the moment 
20 method, this is segmented into elements, and the mutual 
impedance etc. among elements are found by predetermined 
computations . The mutual impedance etc . found and wave 
sources specified by the configuration information (wave 
source of electronic apparatus and wave source of 
25 antenna) are substituted into the simultaneous equations 

of the moment method and the equations are solved so as 
to find the electric current and magnetic current flowing 
through the electronic apparatus. Therefore, it has 
become possible to simulate the effect of a radio wave 
30 radiated from an antenna on an electronic apparatus. 

In view of this, the present inventors 
disclosed in Japanese Patent Application No. 9-90412, 
corresponding to USSN 08/803,166 and German Patent 
Application No. 9710787.0, the use of the moment method 
35 to simulate the effect of a radio wave radiated from an 

antenna on an electronic apparatus . 

In this invention, note was taken of the fact 
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that if the frequency of the carrier wave is f c and the 
frequency of the modulation wave is f m , when amplitude 
modulation is applied, the frequency of the radio wave 
radiated from an antenna can be broken down into three 
parts, that is, f c/ (f c +f ffi ), and (f c -f m ) and the moment 
method is applied by setting the electronic apparatus to 
be simulated and the antenna radiating the radio wave as 
a single object for application of the moment method, 
where the moment method is applied to these above three 
wave sources. Thus it becomes possible to simulate the 
effect of the radio wave radiated from an antenna on the 
electronic apparatus . 

It is true that it is possible to simulate the 
effects of a radio wave radiated from an antenna on an 
electronic apparatus according to the invention of 
Japanese Patent Application No. 9-90412. 

According to Japanese Patent Application No. 9- 
90412, however, it is necessary to calculate the mutual 
impedance, mutual admittance, and mutual reaction and 
solve simultaneous equations under the moment method for 
the frequency f c , to calculate the mutual impedance, 
mutual admittance, and mutual reaction and solve 
simultaneous equations under the moment method for the 
frequency (f c +f m ), and to calculate the mutual impedance, 
mutual admittance, and mutual reaction and solve 
simultaneous equations under the moment method for the 
frequency (f c -f m ). 

This calculation of the mutual impedance, 
mutual admittance, and mutual reaction, however, takes an 
extremely long time. Due to this, there is the problem 
that high speed simulation of the effect of a radio wave 
radiated from an antenna on an electronic apparatus is 
not possible using the invention of Japanese Patent 
Application No. 9-90412. 

Note that in the invention of Japanese Patent 
Application No. 9-90412, the present inventors mainly 



disclosed technology for simulation of an electric 
current, magnetic current, and intensity of an 
electromagnetic field in the time domain using the method 
of high speed calculation of the mutual impedance, mutual 
admittance, and mutual reaction disclosed in Japanese 
Unexamined Patent Publication (Kokai) No. 9-196986 
(Japanese Patent Application No. 7-298062). 

Further, to simulate the effect of a radio wave 
radiated from an antenna on an electronic apparatus, it 
is necessary that the intensity of the electric field 
applied to the electronic apparatus satisfy local legal 
requirements . 

Japanese Patent Application No. 9-90412 does 
not consider this point, however. When using the 
invention of Japanese Patent Application No. 9-90412 to 
simulate the effect of a radio wave radiated from an 
antenna on an electronic apparatus, it is necessary to 
change the positions of the antenna and electronic 
apparatus on a trial and error basis for the simulation. 
There will therefore be the problem that high speed 
simulation of the effect of a radio wave radiated from an 
antenna on an electronic apparatus will not be possible. 
SUMMARY OF THE INVENTION 

The present invention was made in consideration of 
this situation. An object of the present invention is to 
provide a novel apparatus and method for calculating 
immunity from a radiated electromagnetic field which make 
possible high speed simulation of the electric current 
flowing through an electronic apparatus due to a radio 
wave radiated from an antenna and to provide a storage 
medium storing programs used for the same. 

To attain the above object, the present invention 
divides a radio wave radiated from an antenna into a 
carrier wave, upper sideband wave, and lower sideband 
wave and uses the moment method to simulate the effect of 
the radio wave on an electronic apparatus . It calculates 
the mutual impedance for just one frequency component out 
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of the above three frequency components and uses the 
calculated mutual impedance to solve the simultaneous 
equations under the moment method so as to calculate the 
electric current flowing through the electronic apparatus 
5 and uses that mutual impedance to solve the simultaneous 

equations under the moment method for one frequency among 
them, while ignoring the wave source of the electronic 
apparatus, so as to calculate the electric current of 
that frequency component flowing through the electronic 
10 apparatus and calculates the electric currents of the 

remaining frequency components by a proportional 
operation. By this, it is able to calculate the electric 
current flowing through the electronic apparatus due to a 
radio wave radiated from an antenna. As a result, it is 
15 possible to realize an apparatus for calculating immunity 

from a radiated electromagnetic field which makes 
possible high speed simulation of the electric current 
flowing through an electronic apparatus due to a radio 
wave radiated from an antenna. 
2 0 BRIEF DESCRIPTION OF THE DRAWINGS 

The above object and features of the present 
invention will be more apparent from the following 
description of the preferred embodiments with reference 
to the accompanying drawings, wherein: 
25 Figs. 1A, IB and 1C are views of the basic 

configuration of the present invention; 

Fig. 2 is another view of the basic configuration of 
the present invention; 

Fig. 3 is a view of an embodiment of the present 
30 invention; 

Fig. 4 is a view explaining the test regulations; 

Fig. 5 is a view explaining an antenna model; 

Figs . 6A and 6B are views explaining a radio wave 
radiated from an antenna; 
35 Fig. 7 is a view of the flow of processing of a 

program for generating an antenna model; 

Figs. 8A and 8B are other views of the flow of 



processing by a program for generating an antenna model; 

Fig. 9 is a view of the flow of processing by a 
simulation program; 

Fig. 10 is another view of the flow of processing by 
a simulation program; 

Fig. 11 is a view explaining the method of 
calculation of the mutual impedance; 

Figs. 12A and 12B are other views explaining the 
method of calculation of the mutual impedance; 

Fig. 13 is a view explaining simultaneous equations 
under the moment method; 

Fig. 14 is a view explaining the voltage between 
conductors ; 

Figs. 15A, 15B and 15C are other views explaining 
the voltage between conductors; 

Fig. 16 is a view explaining LDU decomposition; 

Fig. 17 is a view explaining LU decomposition; and 

Fig. 18 is a view explaining simultaneous equations 
under the moment method. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Preferred embodiments of the present invention will 
be described next with reference to the related figures. 

Figures 1A, IB, and 1C and Fig. 2 show the basic 
configuration of the present invention. 

In the figure, 1 is an apparatus for calculating 
immunity from a radiated electromagnetic field to which 
the present invention is applied. The apparatus simulates 
the electric current flowing through an electronic 
apparatus due to a radio wave radiated from an antenna. 

The apparatus 1 for calculating immunity from a 
radiated electromagnetic field reads the configuration 
information of the electronic apparatus to be simulated 
from an electronic apparatus configuration data file 2 
and segments the electronic apparatus and the antenna 
used for the simulation (whose configuration information 
is read from an antenna configuration data file 4 shown 
in Fig. 2) into elements. Explaining the case where no 
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consideration is given to a dielectric , the mutual 
impedance among elements is calculated and simultaneous 
equations under the moment method defining the 
relationship among the mutual impedance, wave sources, 
5 and electric currents flowing through elements are solved 

so as to simulate the electric current flowing through an 
electronic apparatus due to a radio wave radiated by an 
antenna. The results of that simulation are output to an 
output device 3 . 
10 The apparatus 1 for calculating immunity from a 

radiated electromagnetic field of the present invention 
of the basic configuration shown in Fig. 1A is provided 

O with a first calculating means 10, a decomposing means 

jfj 11, and a second calculating means 12. 

?X= 15 The first calculating means 10 sets a representative 

frequency (for example, the carrier wave) with respect to 
fil the carrier wave frequency, upper sideband frequency and 

m lower sideband frequency of the radio wave radiated by an 

i«i antenna and calculates the mutual impedance among 

U 20 elements at that representative frequency. 

\"* The decomposing means 11 applies LU decomposition or 

i|3 LDU decomposition to the matrix of the mutual impedance 

calculated by the first calculating means 10. 

The second calculating means 12 solves the 
25 simultaneous equations under the moment method having the 

mutual impedance calculated by the first calculating 
means 10 for the carrier wave frequency, upper sideband 
frequency and lower sideband frequency so as to calculate 
the electric current flowing through the electronic 
30 apparatus due to a radio wave radiated by an antenna. 

When considering a dielectric, the mutual admittance 
and mutual reaction among elements at the representative 
frequency are calculated in addition to the mutual 
impedance and simultaneous equations under the moment 
35 method considering a dielectric having the mutual 

impedance, mutual admittance and mutual reaction are 
solved . 
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The functions of the apparatus 1 for calculating 
immunity from a radiated electromagnetic field of the 
present invention of the basic configuration illustrated 
in Fig. 1A are more specifically realized by programs. 
5 The programs are stored in a floppy disk or other medium, 

stored in the disk etc. of a server etc., and installed, 
from these disks etc., in the apparatus 1 for calculating 
immunity from a radiated electromagnetic field and 
operated in a memory for realization of the present 

10 invention. 

In the apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 
of the basic configuration illustrated in Fig. 1A, the 
first calculating means 10 sets a representative 

15 frequency taking into consideration the fact that there 

is only a slight difference among the carrier wave 
frequency, upper sideband frequency and lower sideband 
frequency forming a radio wave radiated by an antenna. It 
then calculates the mutual impedance among elements at 

2 0 that representative frequency to calculate the mutual 

impedance common to these frequencies. 

Receiving this mutual impedance, the second 
calculating means 12 solves the simultaneous equations 
under the moment method having the calculated mutual 

25 impedance for the carrier wave frequency, upper sideband 

frequency and lower sideband frequency to calculate the 
electric current flowing through the electronic apparatus 
due to a radio wave radiated by an antenna. 

At this time, when the decomposing means 11 is 

30 provided, the second calculating means 12 uses the LU 
decomposed or LDU decomposed matrix of the mutual 
impedance to solve the moment method. LU decomposition or 
LDU decomposition take time, but it is possible to solve 
the moment method at a high speed by using an LU 

35 decomposed or LDU decomposed matrix of the mutual 

impedance. In total, therefore, it becomes possible to 
solve simultaneous equations under the moment method at a 



high speed. 

In this way, in the apparatus 1 for calculating 
immunity from a radiated electromagnetic field of the 
present invention of the basic configuration shown in 
Fig. 1A, when configured to break down a radio wave 
radiated by an antenna into a carrier wave, upper 
sideband wave and lower sideband wave and use the moment 
method to simulate the effect of the radio wave radiated 
by an antenna, the mutual impedance is calculated for 
just one frequency component among the carrier wave 
frequency, upper sideband frequency and lower sideband 
frequency, the calculated mutual impedance is used to 
solve the simultaneous equations under the moment method 
for the carrier wave frequency so as to calculate the 
electric current of the carrier wave frequency component 
flowing through the electronic apparatus due to a radio 
wave radiated by an antenna, to solve the simultaneous 
equations under the moment method for the upper sideband 
frequency so as to calculate the electric current of the 
upper sideband frequency component flowing through the 
electronic apparatus due to a radio wave radiated by an 
antenna, and to solve the simultaneous equations under 
the moment method for the lower sideband frequency so as 
to calculate the electric current of the lower sideband 
frequency component flowing through the electronic 
apparatus due to a radio wave radiated by an antenna. 
Therefore, it is possible to simulate the current flowing 
through the electric current due to a radio wave radiated 
by an antenna at a high speed. 

The apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 
of the basic configuration shown in Fig. IB is provided 
with a first calculating means 20, a second calculating 
means 21, and a third calculating means 22. 

The first calculating means 20 sets a representative 
frequency out of the carrier wave frequency, upper 
sideband frequency, and lower sideband frequency forming 
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a radio wave radiated by an antenna and calculates the 
mutual impedance among elements at that representative 
frequency. 

The second calculating means 21 solves the 
5 simultaneous equations under the moment method having the 

mutual impedance calculated by the first calculating 
means 20 for one of the carrier wave frequency, upper 
sideband frequency and lower sideband frequency, while 
ignoring the wave source of the electronic apparatus, so 
10 as to calculate the electric current flowing through the 

electronic apparatus due to a radio wave radiated from an 
antenna . 

O The third calculating means 22 calculates the 

ifi electric currents, other than the electric current 

n:l 15 calculated by the second calculating means 21, flowing 

through the electronic apparatus due to a radio wave 
Hj radiated from an antenna, by a proportional operation, by 

Ui using the electric current calculated by the second 

p calculating means 21 and the value of the wave source of 

2 0 the antenna . 

jT When considering a dielectric, the mutual admittance 

*S and mutual reaction among elements at the representative 

frequency are calculated in addition to the mutual 

impedance and simultaneous equations under the moment 
25 method, considering a dielectric having the mutual 

impedance, mutual admittance and mutual reaction, are 

solved . 

The functions of the apparatus 1 for calculating 
immunity from a radiated electromagnetic field of the 

30 present invention of the basic configuration illustrated 

in Fig. IB are more specifically realized by programs. 
The programs are stored in a floppy disk or other medium, 
stored in the disk etc. of a server etc., and installed, 
from these disks etc . , in the apparatus for calculating 

35 immunity from a radiated electromagnetic field 1 and 

operated in a memory for realization of the present 
invention . 
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In the apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 
of the basic configuration illustrated in Fig. IB, the 
first calculating means 20 sets a representative 
5 frequency taking into consideration the fact that there 

is only a slight difference among the carrier wave 
frequency, upper sideband frequency and lower sideband 
frequency forming a radio wave radiated by an antenna. It 
then calculates the mutual impedance among elements at 
10 that representative frequency to calculate the mutual 

impedance common to these frequencies. 

Receiving this mutual impedance, explaining the case 
where the second calculating means 21 uses the carrier 
flj wave frequency in its calculations, the second 

15 calculating means 21 solves the simultaneous equations 

J:; under the moment method for the carrier wave frequency 

jy having the calculated mutual impedance, while ignoring 

./ the wave source of the electronic apparatus, to calculate 

O the electric current of the carrier wave frequency 

J") 2 0 component flowing through the electronic apparatus due to 

Hi a radio wave radiated by an antenna. 

Receiving the electric current of the carrier wave 
frequency component calculated by the second calculating 
means 21, the third calculating means 22 calculates the 
2 5 electric current of the upper sideband frequency 

component flowing through the electronic apparatus due to 
a radio wave radiated by an antenna, by a proportional 
operation, by using the calculated electric current of 
the carrier wave frequency component, the value of the 
30 wave source of the antenna at the carrier wave frequency 
and the value of the wave source of the antenna at the 
upper sideband frequency and also calculates the electric 
current of the lower sideband frequency component flowing 
through the electronic apparatus due to a radio wave 
35 radiated by an antenna, by a proportional operation, by 

using the calculated electric current of the carrier wave 
frequency component, the value of the wave source of the 
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antenna at the carrier wave frequency and the value of 
the wave source of the antenna at the lower sideband 
frequency. 

In this way, in the apparatus 1 for calculating 
immunity from a radiated electromagnetic field of the 
present invention of the basic configuration shown in 
Fig. IB, when configured to break down a radio wave 
radiated by an antenna into a carrier wave, upper 
sideband wave and lower sideband wave and use the moment 
method to simulate the effect of a radio wave radiated by 
an antenna, the mutual impedance is calculated for just 
one frequency component among the above three 
frequencies, the calculated mutual impedance is used to 
solve the simultaneous equations under the moment method 
for one of the frequencies among these, while ignoring 
the wave source of the electronic apparatus, so as to 
calculate the electric current of that frequency 
component flowing through the electronic apparatus, and 
the electric currents of the remaining frequency 
components are calculated by a proportional operation, so 
it is possible to simulate the current flowing through 
the electric current due to a radio wave radiated by an 
antenna at a high speed. 

The apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 
of the basic configuration shown in Fig. 1C is provided 
with a first calculating means 30, a decomposing means 
31, a second calculating means 32, a third calculating 
means 33, and a fourth calculating means 34. 

The first calculating means 30 sets a representative 
frequency out of the carrier wave frequency, upper 
sideband frequency and lower sideband frequency forming a 
radio wave radiated by an antenna and calculates the 
mutual impedance among elements at that representative 
frequency. 

The decomposing means 31 applies LU decomposition or 
LDU decomposition on the matrix of the mutual impedance 
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calculated by the first calculating means 30. 

The second calculating means 32 solves the 
simultaneous equations under the moment method having the 
mutual impedance calculated by the first calculating 
5 means 30 for the frequency, among the carrier wave 

frequency, upper sideband frequency and lower sideband 
frequency which overlaps the frequency, including a 
higher harmonic component, of the wave source of the 
electronic apparatus, so as to calculate the electric 

10 current flowing through the electronic apparatus due to a 

radio wave radiated by an antenna. 

The third calculating means 33 solves the 
simultaneous equations under the moment method having the 
mutual impedance calculated by the first calculating 

15 means 30 for one of the frequencies not used in the 

calculation by the second calculating means 32 so as to 
calculate the electric current, other than the electric 
current calculated by the second calculating means 32, 
flowing through the electronic apparatus due to a radio 

20 wave radiated from an antenna. 

The fourth calculating means 34 calculates the 
electric current, other than the electric currents 
calculated by the second and third calculating means 32 
and 33, flowing through an electronic apparatus due to a 

25 radio wave radiated from an antenna, by a proportional 

operation, by using the electric current calculated by 
the third calculating means 33 and the value of the wave 
source of the antenna. 

When considering a dielectric, the mutual admittance 

30 and mutual reaction among elements at the representative 

frequency are calculated in addition to the mutual 
impedance and simultaneous equations under the moment 
method having the mutual impedance, mutual admittance and 
mutual reaction, considering a dielectric, are solved. 

35 The functions of the apparatus 1 for calculating 

immunity from a radiated electromagnetic field of the 
present invention of the basic configuration illustrated 
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in Fig. 1C are more specifically realized by programs. 
The programs are stored in a floppy disk or other medium, 
stored in the disk etc. of a server etc., and installed, 
from these disks etc . , in the apparatus 1 for calculating 
5 immunity from a radiated electromagnetic field and 

operated in a memory for realization of the present 
invention . 

In the apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 

10 of the basic configuration illustrated in Fig. 1C, the 

first calculating means 30 sets a representative 
frequency taking into consideration the fact that there 
is only a slight difference among the carrier wave 
frequency, upper sideband frequency and lower sideband 

15 frequency forming a radio wave radiated by an antenna. It 

then calculates the mutual impedance among elements at 
that representative frequency to calculate the mutual 
impedance common to these frequencies . 

Receiving this mutual impedance, explaining the case 

20 where the carrier wave frequency overlaps the frequency 
of the wave source of the electronic apparatus and the 
upper sideband frequency and lower sideband frequency do 
not overlap the frequency of the wave source of the 
electronic apparatus, the second calculating means 32 

2 5 solves the simultaneous equations under the moment method 

having the calculated mutual impedance (wave source of 
electronic apparatus is considered), for the carrier wave 
frequency, to calculate the electric current of the 
carrier wave frequency component flowing through the 

30 electronic apparatus due to a radio wave radiated by an 

antenna, while the third calculating means 33 solves the 
simultaneous equations under the moment method having the 
calculated mutual impedance (wave source of electronic 
apparatus is not considered) , for example for the upper 

35 sideband frequency, to calculate the electric current of 

for example the upper sideband frequency component 
flowing through the electronic apparatus due to a radio 
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wave radiated by an antenna. 

Receiving the electric current of the upper sideband 
frequency component calculated by the third calculating 
means 3 3 , the fourth calculating means 34 calculates the 
5 electric current of the lower sideband frequency 

component flowing through the electronic apparatus due to 
a radio wave radiated from an antenna, by a proportional 
operation, by using the calculated electric current of 
the upper sideband frequency component, the value of the 
10 wave source of the antenna at the upper sideband 

frequency and the value of the wave source of the antenna 
at the lower sideband frequency. 
O At this time, when the decomposing means 31 is 

P provided, the second and third calculating means 32 and 

(33 15 33 solve the simultaneous equations of the moment method 

v i using the LU decomposed or LDU decomposed matrix of the 

Hj mutual impedance. LU decomposition and LDU decomposition 

— take time, but it is possible to solve the simultaneous 

ri equations under the moment method using a LU decomposed 

l*J 2 0 or LDU decomposed matrix of mutual impedance at a high 

Yl speed. In total, therefore, it becomes possible to solve 

>Q the simultaneous equations under the moment method at a 

it 1 : 

high speed. 

In this way, in the apparatus 1 for calculating 
2 5 immunity from a radiated electromagnetic field of the 

present invention of the basic configuration shown in 
Fig. 1C, when configured to break down a radio wave 
radiated by an antenna into a carrier wave, upper 
sideband wave and lower sideband wave and use the moment 
30 method to simulate the effect of a radio wave radiated by 

an antenna, the mutual impedance is calculated for just 
one frequency component among the above three 
frequencies, the calculated mutual impedance is used to 
solve the simultaneous equations under the moment method 
35 for the frequency overlapping the frequency of the wave 

source of the electronic apparatus so as to calculate the 
electric current of that frequency component flowing 
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through the electronic apparatus, the simultaneous 
equations under the moment method are solved for just one 
of the electric currents of the non-overlapping frequency- 
components to calculate the electric current of that 
5 frequency component flowing through the electronic 

apparatus, and the electric current of the remaining 
frequency component is calculated by a proportional 
operation, so it is possible to simulate the current 
flowing through the electric current due to the radio 

10 wave radiated by an antenna at a high speed. 

The apparatus 1 for calculating immunity from a 
radiated electromagnetic field of the present invention 
of the basic configuration shown in Fig. 2 is provided 
with a managing means 40, a first computing means 41, a 

15 second computing means 42, an executing means 43, an 

acquiring means 44, a calculating means 45, a setting 
means 46, and an alarm means 47. 

The managing means 40 manages antenna information 
realizing a prescribed intensity of an electric field on 

20 an electronic apparatus (intensity of electric field 

prescribed in regulations etc . ) 

The first computing means 41 assumes a state where 
there is no electronic apparatus, segments the antenna to 
be registered in the managing means 40 into elements, 

25 calculates the mutual impedance among these elements, and 

solves the simultaneous equations under the moment method 
defining the relationship among the calculated mutual 
impedance, wave source of the antenna and electric 
currents flowing through the elements so as to calculate 

30 the electric currents flowing through these antenna 

elements. In this calculation processing, the first 
computing means 41 solves the simultaneous equations 
under the moment method for one frequency among the 
carrier wave frequency, upper sideband frequency and 

35 lower sideband frequency so as to calculate the electric 

current flowing through the antenna. 

The second computing means 42 calculates the 
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intensity of the electric field which the electric 
current calculated by the first computing means 41 causes 
in the electronic apparatus at different locations of 
installation . 

5 The executing means 4 3 changes the distance between 

the antenna and electronic apparatus and the value of the 
wave source of the antenna to determine the specific 
distance and value of the wave source giving the 
prescribed intensity of electric field calculated by the 
10 second computing means 42 and registers the thus 

prescribed antenna information in the managing means 40. 

The acquiring means 44 acquires an antenna model for 
□ use in the simulation from the managing means 40 when a 

:H simulation request is issued. 

!jj 15 The calculating means 45 segments the electronic 

apparatus and the antenna designated by the antenna 
;y information acquired by the acquiring means 44 into 

ul elements, calculates the mutual impedance among these 

■«! elements, and solves the simultaneous equations under the 

Ijj 20 moment method defining the relationship among the 

Y": calculated mutual impedance, wave source and electric 

currents flowing through the elements so as to calculate 
the electric current flowing through the electronic 
apparatus due to a radio wave radiated by an antenna. In 
25 this calculation processing, the calculating means 45 can 

calculate the electric current using the high speed 
calculation techniques executed by the apparatuses 1 for 
calculating immunity from a radiated electromagnetic 
field of the present invention of the basic 
30 configurations shown in Figs. 1A to 1C . 

The setting means 46 sets a threshold current for a 
specified element and sets a threshold voltage for a 
position between specified conductor elements. 

The alarm means 4 7 outputs information on whether 
35 the electric current flowing through a specified element 

exceeds a threshold current set by the setting means 46 
and outputs information on whether the voltage generated 
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at a position between specified conductor elements 
exceeds a threshold voltage set by the setting means 46. 

When considering a dielectric, the mutual admittance 
and mutual reaction among elements are calculated in 
5 addition to the mutual impedance and simultaneous 

equations under the moment method having the mutual 
impedance, mutual admittance and mutual reaction, 
considering a dielectric, are solved. 

The functions of the apparatus 1 for calculating 

10 immunity from a radiated electromagnetic field of the 

present invention of the basic configuration illustrated 
in Fig. 2 are more specifically realized by programs. The 
programs are stored in a floppy disk or other medium, 
stored in the disk etc. of a server etc., and installed, 

15 from these disks etc., in the apparatus 1 for calculating 

immunity from a radiated electromagnetic field and 
operated in a memory for realization of the present 
invention . 

In the apparatus 1 for calculating immunity from a 
20 radiated electromagnetic field of the present invention 

of the basic configuration illustrated in Fig. 2, when 
the first computing means 41 assumes a state where there 
is no electronic apparatus to calculate the electric 
currents flowing through elements of the antenna, the 
25 second computing means 42 calculates the intensity of the 

electric field caused by the calculated electric current 
at the electric current at different locations of 
installation . 

Receiving the calculation processing of the second 
30 computing means 42, the executing means 43 changes the 

distance between the antenna and electronic apparatus and 
the value of the wave source of the antenna to determine 
the specific distance and value of the wave source giving 
the prescribed intensity of electric field calculated by 
35 the second computing means 42 and registers the thus 

prescribed antenna information in the managing means 40. 
Receiving the antenna information managed by the 
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managing means 40, the acquiring means 44 acquires 
antenna information for use in the simulation from the 
managing means 40 when a simulation request is issued. 
Receiving this, the calculating means 45 segments the 
5 electronic apparatus to be simulated and the antenna 

designated by the antenna information (Fig. 5) acquired 
by the acquiring means 44 into elements, calculates the 
mutual impedance among these elements, and solves the 
simultaneous equations under the moment method defining 

10 the relationship among the calculated mutual impedance, 

wave source and electric currents flowing through the 
elements so as to calculate the electric current flowing 
through the electronic apparatus due to a radio wave 
radiated by an antenna. 

15 At this time, the alarm means 47 compares the 

electric current, flowing through a specified element, 
calculated by the calculating means 45 and the threshold 
current set by the setting means 4 6 and outputs 
information on whether the electric current exceeds the 

2 0 threshold current. Further, the alarm means 4 7 compares 

(i) the voltage generated at a position between specified 
conductor elements, derived by making the voltage 
generated across a resistor, virtually inserted between 
the conductors, one obtained if the resistor has an 

25 infinitely large resistance, and (ii) the threshold 

voltage set by the setting means 4 6 and outputs 
information on whether this voltage exceeds the threshold 
voltage . 

In this way, in the apparatus 1 for calculating 
30 immunity from a radiated electromagnetic field of the 

present invention of the basic configuration shown in 
Fig. 2, when configured to use the moment method to 
simulate the effect of a radio wave radiated by an 
antenna, antenna information for realizing a prescribed 
35 intensity of an electric field on the electronic 

apparatus (intensity of electric field prescribed in 
regulations etc . ) is prepared in advance and the antenna 
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information is used to solve the simultaneous equations 
under the moment method so as to calculate the electric 
current flowing through the electronic apparatus when 
there is a simulation request, so it is possible to 
5 simulate at a high speed the electric current flowing 

through the electronic apparatus when a prescribed 
intensity of electric field is applied due to a radio 
wave radiated by an antenna . 

The present invention will be explained in further 
10 detail below in accordance with specific embodiments. 

Figure 3 is a view of an embodiment of the apparatus 
1 for calculating immunity from a radiated 
O electromagnetic field of the present invention. 

The apparatus 1 for calculating immunity from a 
m 15 radiated electromagnetic field of the present invention 

^i; according to this embodiment is provided with an 

ilj electronic apparatus configuration data file 2 for 

* ! managing information on the configurations of electronic 

p apparatuses to be simulated, an output device 3 for 

1=4 20 outputting the results of the simulation, an antenna 
^2 configuration data file 4 for managing information on the 

=D configurations of antennas used for the processing in the 

*' s " simulation, an antenna model generation program 100, 

installed through a floppy disk or through a line etc., 
25 for reading the configuration information of an antenna 

from the antenna configuration information data file 4 
and generating an antenna model used for the processing 
in the simulation, an antenna model library 200 for 
managing the antenna models generated by the antenna 
30 model generation program 100, and a simulation program 
300, installed through a floppy disk or through a line 
etc., for reading the configuration information of an 
electronic apparatus to be simulated from the electronic 
apparatus configuration data file 2, simulating the 
35 electric current flowing through the electronic apparatus 

due to a radio wave radiated by an antenna, and 
outputting the result of the simulation to the output 
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device 3 . 

A test for examining the effect of a radio wave 
radiated by an antenna on an electronic apparatus is 
performed, as shown in Fig. 4, by establishing a testing 
5 zone having a size of for example 1.5 m x 1.5 m in which 

an electronic apparatus is arranged at a position of a 
certain height from the ground (set in a state 
perpendicular to the ground as shown in the figure) and 
applying an electric field of a prescribed magnitude to 
10 the testing zone by a radio wave radiated by an antenna 

(for example, an electric field of 3V/m with a difference 
between the maximum value and minimum value of not more 
^ than 6 dB) . 

The antenna model generation program 100 generates 
is an antenna model realizing the test conditions. 

Specifically, as shown in Fig. 5, it generates an antenna 
jlj model comprised of a configuration ID (specifying the 

yi antenna configuration information stored in the antenna 

□ configuration data file 4), the class of the antenna 

j«? 20 (dipole, Log-Peri, biconical, etc.), the distance between 

j"]: the antenna and the electronic apparatus, the height of 

?S the antenna, the modulation condition of the antenna, the 

' hi direction in which the electric field is to be applied 

from the antenna (horizontal direction/vertical 
25 direction, front/rear/right/left), the allowable level of 

the intensity of the electric field in the testing zone 
(above-mentioned 6 dB, 3V/m) , and the location of the 
uniform electric field plane of the testing zone (above- 
mentioned 1.5 m x 1.5 m) . 
30 Here, the modulation conditions of the antenna 

specifically consist of information (condition) on the 
frequency range of the carrier wave, the frequency f c of 
the carrier wave, the frequency f m of the modulation 
wave, and whether the modulation is amplitude modulation 
35 or pulse modulation. 

The radio wave radiated by an antenna, if the 



modulation signal is v(t) f becomes: 

f(t) = Ao[l+kv(t) ]exp( j(o c t) 
The modulation signal v(t) is expanded by a complex 
Fourier series as follows: 

v(t) = Ec n exp( jnw m t) n=0±l±2... 

In this equation, in the case of amplitude 
modulation, "n=l". On the other hand, in the case of 
pulse modulation, "n=0 to +«>, but in the simulation 
processing, "n=0 to ±L" is set. Therefore, 

v(t) = £c n exp( jno) m t) n=0±l...±L 
If this is substituted in 

f(t) = Ao[l+kv(t) ]exp( ja) c t) 
the result is , 

f(t) = Aq[ jo c t)+AokEc n exp[ j(o c +nw m )t] 

From this, in the case of amplitude modulation, as 
shown in Fig. 6A, the wave source of the antenna can be 
broken down into three waves sources, that is, the wave 
source of the frequency f c of the carrier wave, the wave 
source of the frequency (f c +f m ) of the upper sideband 
wave and the wave source of the frequency (f c -f m ) of the 
lower sideband wave. Further, in the case of pulse 
modulation, as shown in Fig. 6B, it may be broken down 
into the three types of wave sources of the wave source 
of the frequency f c of the carrier wave, the wave sources 
of the frequencies increased incrementally by f m from f c/ 
which become the frequency of the upper sideband wave, 
and the wave sources of the frequencies decreased 
incrementally by f m from f c , which become the frequency 
of the lower sideband wave. 

From this, in the antenna model, information on the 
frequency range of the carrier wave, the frequency f c of 
the carrier wave, the frequency f m of the modulation 
wave, and whether the modulation is amplitude modulation 
or pulse modulation is managed. 

Figure 7 and Figs . 8A and 8B show an embodiment of 
the flow of processing executed by the antenna model 
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generation program 100, while Fig. 9 and Fig. 10 show an 
example of the flow of processing executed by the 
simulation program 300. Next, the present invention will 
be explained in detail in accordance with these flows of 
5 processing . 

First, an explanation will be made of the processing 
executed by the antenna model generation program in 
accordance with the flow of processing of Fig. 7 and 
Figs . 8A and 8B. 

10 When an operator specifies a configuration ID and 

issues a request for generation of an antenna model, the 
antenna model generation program 100, as shown in the 
flow of processing of Fig. 7 and Figs. 8A and 8B, first, 
at step (ST) 1, receives as input the configuration 

15 information of the antenna indicated by the specified 

configuration ID from the antenna configuration data file 
4, then, at step 2, segments the input antenna to which 
the simultaneous equations under the moment method are to 
be applied, into elements. 

20 Next, at step 3, it receives as input the test 

conditions for the simulation for the antenna through a 
dialog with the operator. That is, it receives as input 
the frequency range of the carrier wave from the antenna 
(for example, 30 MHz to 1 GHz), the frequency of the 

25 modulation wave from the antenna (for example, 1 kHz), 

the modulation mode (amplitude modulation/pulse 
modulation), the height of the antenna, the direction to 
be radiated from the antenna, the size of the testing 
zone (for example, 1.5 m x 1.5 m), the allowable level of 

30 the intensity of the electric field in the testing zone 

(for example, 6 dB and 3 V/m) , and other information 
constituting the antenna model . 

Here, the distance between the antenna and the 
testing zone is determined from now so as to realize the 

35 input allowable level of the intensity of the wave source 

of the antenna, so the prescribed initial value is input. 
At this time, looking at the initial value of the 
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distance between the antenna and testing zone, as will be 
understood from the later explanation, a small value is 
set for the distance used in the actual simulation as 
well . 

Next, at step 4, the program selects one carrier 
wave frequency in a manner increased incrementally for 
example by 50 MHz from the frequency range of the carrier 
wave input at step 3. Next, at step 5, it judges if all 
of the frequencies have finished being selected. When it 
judges that not all of the frequencies have finished 
being selected, that is, when it judges that a carrier 
wave frequency could be selected at step 4, it proceeds 
to step 6 , where it calculates the mutual impedance among 
elements segmented at step 2 at the selected carrier wave 
frequency. 

The mutual impedance is specifically calculated by 
assuming the monopoles as shown in Fig. 11 (<1> to <4> in 
the figure) . 

The general equation of the mutual impedance 
between an element i and an element j is shown by the 
equation illustrated in Fig. 12A. In the figure, a> 
indicates the angular frequency, r the distance between 
monopoles, Ji and J 2 the shapes of distribution of 
electric current on the monopoles , and <t> the angle 
between monopoles and pi=( -1/ jo)x[3J 1 /5t ] and p 2 =(- 
1/ j(D)x[dJ 2 /dt] . 

As the distributions J lr J 2 of electric current on 
the monopoles, assuming that: 

Electric current monopole <1> J x = sink( z-z 0 ) /sinkdi 
Electric current monopole <2> J x = sink{ z 2 -z ) /sinkd 2 
Electric current monopole <3> J 2 = sink( t-t 0 ) /sinkd 3 
Electric current monopole <4> J 2 = sink( t 2 -t ) /sinkd^ 
where, d x is the length of the monopole <1>, d 2 is 
the length of the monopole <2> , d 3 is the length of the 
monopole <3>, and d A is the length of the monopole <4>, 
the mutual impedance Z 13 between the monopole <1> 
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and the monopole <3> and the mutual impedance Z 14 between 
the monopole <1> and the monopole <4> become as given by 
the equations illustrated in Fig. 12B. The mutual 
impedance Z 23 between the monopole <2> and the monopole 
5 <3> and the mutual impedance Z 24 between the monopole <2> 
and the monopole <4> are given by similar equations. 

Next, at step 6 (Fig. 7), the program computes these 
equations to find the mutual impedance Z ±j ( =Z 13 +Z i4+Z 23 +Z 2/( ) 
between the element i and the element j . 
10 When calculating at step 6 the mutual impedance 

between elements at the carrier wave frequency selected 
at step 4, the program then solves, at step 7, the 
;ji simultaneous equations under the moment method, shown by 

jlj the equation in Fig. 13, having that mutual impedance so 

15 as to calculate the electric current flowing through the 

I 'J 

Ji; antenna elements segmented at step 2 . 

Next, at step 8, the program computes the intensity 
of the electric field applied by the electric current 
□ calculated at step 7 to a plurality of measurement points 

J 5 /} 2 0 (for example, 16 points) set in the testing zone and 

!»a input at step 3. The processing for this calculation is 

executed in accordance with the known Maxwell's 
electromagnetic theoretical equations. Next, at step 9, 
the program specifies the maximum value and minimum value 
25 of the intensity of the electric field calculated at step 

8 and calculates the difference between these two 
intensities of electric field to check for the degree of 
uniformity of the electric field in the testing zone. 
Next, at step 10, the program checks if the 
30 difference of the intensities of the electric field 

calculated at step 9 is less than an allowable 
attenuation ratio (for example, 6 dB) input at step 3. 
When it is judged by this checking processing that it is 
not less than the allowable attenuation ratio, the 
35 program proceeds to step 11, where it increases the 

distance between the antenna and testing zone so that the 
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difference becomes less than the allowable attenuation 
ratio. Note that when increasing the distance, the 
uniformity of the electric field in the testing zone is 
improved, so it is not necessary to repeat the checking 
5 processing for the carrier wave frequencies for which the 

check of uniformity of the electric field had been 
finished up to then. 

When judging at step 10 that the difference is less 
than the allowable attenuation ratio or when judging at 

10 step 10 that the difference is not less than the 

allowable attenuation ratio and changing the distance at 
step 11 so that the distance becomes less than the 
allowable attenuation ratio, the program proceeds to step 
12 (flow of processing of Fig. 8) where it compares the 

15 minimum value of the intensity of the electric field 

calculated presently at step 8 and the minimum value of 
the intensity of the electric field calculated in the 
past (held in the working area (see WA in Fig. 3)). 

Next, at step 13, the program judges if the minimum 

20 value of the intensity of the electric field currently 

calculated is smaller than that of the past by the 
comparison processing at step 12. If judging that it is 
smaller, the program proceeds to step 14, where it 
updates the minimum value of the intensity of the 

25 electric field held at the above working area, then 

returns to step 4 for the processing of the next carrier 
wave frequency. Conversely, if judging that it is larger, 
the program returns to step 4 for the processing of the 
next carrier wave frequency without performing the 

30 processing of step 14. 

By repeating the processings of step 4 to step 14, 
the distance between the antenna and the testing zone 
giving a difference within the allowable attenuation 
ratio (for example, 6 dB) is found. When the distance 

35 finishes being determined by judging at step 5 that all 

of the frequencies have finished being selected, the 
program proceeds to step 15, where it determines the 
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value of the wave source of the antenna to be used for 
the simulation by using the minimum value of the 
intensity of the electric field held in the working area, 
the value of the wave source of the antenna assumed as 
5 the initial value and the intensity of the electric field 

of the test conditions input at step 3 (for example, 
3V/m) . 

For example, when the intensity of the electric 
field of the test conditions input at step 3 is 3 V/m and 
10 the program assumes the value of the wave source of the 

antenna to be IV (originally expressed as a complex 
number since a phase is involved) to determine the 
distance between the antenna and the testing zone, when 
the minimum value of the intensity of the electric field 
15 held in the working area is for example 1.5 V/m, the 

program determines the value of the wave source of the 
antenna to be used for the simulation to be 2V in 
accordance with the proportional operation of: 
g IV x [3+1.5] = 2V 

kj 2 0 Further, finally, at step 16, the program generates 

an antenna model in accordance with the information of 

=n the test conditions input at step 3, the distance between 

the antenna and testing zone finally determined in 
accordance with the processing of step 11 and the value 
2 5 of the wave source determined at step 15, registers this 

in the antenna model library 200, and then ends the 
processing. 

In this way, in accordance with the flow of 
processing of Fig. 7 and Figs. 8A and 8B, the antenna 

30 model generation program 100 generates an antenna model 
having the data configuration shown in Fig. 5 and 
registers this in the antenna model library 200. 

In the flow of processing, the antenna model 
generation program 100 solves the simultaneous equations 

35 under the moment method for the carrier wave frequency so 

as to evaluate the intensity of the electric field 
applied by the electric current of that frequency 
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component flowing through the antenna and thereby 
determine the distance between the antenna and testing 
zone, but it is also possible to solve the simultaneous 
equations under the moment method for the upper sideband 
5 frequency or lower sideband frequency so as to evaluate 

the intensity of the electric field applied by the 
electric current of those frequency components flowing 
through the antenna and determine the distance between 
the antenna and testing zone. 

10 Next, an explanation will be made of the processing 

executed by the simulation program 300 in accordance with 
the flow of processing of Fig. 9 and Fig. 10. Here, in 
the flow of processing, for convenience in the 
explanation, amplitude modulation is used. Accordingly, 

15 as shown in Fig. 6A, one upper sideband wave and one 

lower sideband wave are assumed. 

When an operator specifies an electronic apparatus 
for simulation and a configuration ID of an antenna model 
for use in the simulation and issues a request for 

20 simulation, the simulation program 300, as shown by the 

flow of processing of Fig. 9 and Fig. 10, first, at step 
(ST) 1, receives as input the configuration information 
of the electronic apparatus to be simulated from the 
electronic apparatus configuration data file 2 . 

25 Next, at step 2, it receives as input the antenna 

model indicated by the specified configuration ID from 
the antenna model library 200, then, at step 3, receives 
as input the configuration information of the antenna 
indicated by the specified configuration ID from the 

30 antenna configuration data file 4. Next, at step 4, it 

segments the input electronic apparatus and antenna, to 
which the simultaneous equations under the moment method 
are to be applied, into elements. 

Next, at step 5, the program selects one carrier 

35 wave frequency from the frequency range of the carrier 

wave specified in the antenna model in a manner increased 
incrementally for example by 30 MHz. Next, at step 6, it 
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judges if all of the frequencies have finished being 
selected. When it judges that all of the frequencies have 
finished being selected, it proceeds to step 7 , where it 
outputs the results of the simulation to the output 
5 device 3 and ends the processing. 

On the other hand, when it judges at step 6 that not 
all of the frequencies have finished being selected, that 
is, when it judges that a carrier wave frequency could be 
selected at step 5, it proceeds to step 8, where it 
10 calculates the mutual impedance among elements segmented 
at step 4 at the selected carrier wave frequency in 
accordance with the above calculation technique. 
□ Next, at step 9, the program judges if the wave 

~\ source of the electronic apparatus to be simulated need 

83 15 not be considered. When judging that it need not be 
% considered, that is, when judging that the frequency, 

jlj including a higher harmonic component, of the wave source 

^ x of the electronic apparatus does not overlap the selected 

p carrier wave frequency, or that it overlaps it, but the 

y 20 value of the wave source is small and therefore the wave 
Y'l source can be ignored, the program proceeds to step 10, 

=n where it solves the simultaneous equations, for the 

* u carrier wave frequency, under the moment method (only 

wave source of antenna exists) having the mutual 
25 impedance calculated at step 8 so as to calculate the 

electric current of the carrier wave frequency component 
flowing through the elements of the electronic apparatus 
segmented at step 4 . 

Next, the program proceeds to step 11 (flow of 
30 processing of Fig. 10) where it calculates the electric 

current of the upper sideband frequency component flowing 
through the elements of the electronic apparatus 
segmented at step 4 , in accordance with a proportional 
operation, by using the electric current of the carrier 
35 wave frequency component calculated at step 10, the value 
of the wave source of the antenna at the carrier wave 
frequency and the value of the wave source of the antenna 
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at the upper sideband frequency and calculates the 
electric current of the lower sideband frequency 
component flowing through the elements of the electronic 
apparatus segmented at step 4, in accordance with a 
5 proportional operation, by using the electric current of 

the carrier wave frequency component calculated at step 
10, the value of the wave source of the antenna at the 
carrier wave frequency and the value of the wave source 
of the antenna at the lower sideband frequency. 
10 That is, expressing the electric current of the 

carrier wave frequency component flowing through the 
elements of the electronic apparatus calculated at step 
£Jj 10 by I c , the electric current of the upper sideband 

ill frequency component flowing through the elements of the 

!;*-: 15 electronic apparatus by I H/ the electric current of the 

j\: lower sideband frequency component flowing through the 

j"y elements of the electronic apparatus by I L , the value V c 

.7 of the wave source of the antenna at the carrier wave 

p| frequency by "V c =a c +jb c ", the value V H of the wave source 

jjj 20 of the antenna at the upper sideband frequency by "V H = 

^; a H +jb H ", and the value V L of the wave source of the 

Cf\ antenna at the lower sideband frequency by "V L = a L + jb L " , 

the electric current of the upper sideband frequency 
component flowing through the elements of the electronic 
25 apparatus is calculated in accordance with the 

proportional operation as shown below: 

I H = I c x [ (a H +jb H )/(a c +jb c ) ] 
and the electric current of the lower sideband frequency 
component flowing through the elements of the electronic 
30 apparatus is calculated in accordance with the 

proportional operation as shown below: 
I L = l c x [ (a L +jb L )/(a c +jb c ) ] 

Next, at step 15, the program uses the electric 
currents found at step 10 and step 11 to calculate the 
35 voltage generated at a position between the conductor 
elements specified by the operator. 
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The processing for this calculation is executed by 
calculating the voltage V p (g>) between conductor elements 
in accordance with: 

V p (o) = -El^wJZ^tu) 

where E is n = 1 to M 
when the position p between conductor elements is 
specified,, the electric current flowing through an 
element n is indicated by I n (w) , and the mutual impedance 
between the position p, between conductor elements, and 
the element n is indicated by Z pn («) . 

Explaining this calculation formula, as shown in 
Fig. 14 , if a resistor R is inserted between the 
conductor of the conductor element pi and the conductor 
element p2 , then based on the boundary condition that the 
electric field on the conductors becoming zero, the 
equation shown in Fig. 15A stands. From this, the 
electric current I p between conductors is found by the 
equation shown in Fig. 15B and, from this, the voltage V p 
between conductors is found by the equation shown in Fig. 
15C. In actuality, since no current flows between 
conductors, in the equation shown in Fig. 15C, "R - 00 , 
Ipif ! P 2 - 0" . From this the equation (Fig. 15C) is found. 

That is, this equation is derived by making the 
voltage generated across a resistor, virtually inserted 
between the conductors, one obtained if the resistor has 
an infinitely large resistance. 

At step 15, the program calculates the voltage 
generated at a position between specified conductor 
elements. Next, at step 16, it judges if the electric 
current flowing through an element specified by the 
operator exceeds a prescribed threshold and if the 
voltage generated at the position between conductor 
elements calculated at step 15 exceeds a prescribed 
threshold, registers the results of its judgement, then 
returns to step 5 for the processing of the next carrier 
wave frequency . 
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Here, the thresholds are set for example through a 
dialog with the operator such that, if the thresholds are 
exceeded, there is a possibility of an electronic circuit 
component positioned at the specified element 
5 malfunctioning due to noise. 

On the other hand, when the program judges at step 9 
that the wave source of the electronic apparatus being 
simulated must be considered, the program proceeds to 
step 13, where it applies LDU decomposition on the mutual 
10 impedance Zfz^) calculated at step 8 in accordance with 

the LDU decomposition rules for a matrix. That is, it 
applies LDU decomposition to the mutual impedance Z(z^) 
as shown in Fig. 16. Here, the following stand for the 
i?j matrix D(d ij ) and the matrix L(l ij ): 

K 15 d Li = z iL - Ed kk l ik 2 

i y 

where, E is the sum for k = 1 to (i-1), i = 1 to n 
lij = [z±j - Ed kk li k ljk]/djj 

where, E is the sum for k = 1 to (j-1), i = 1 ton, 

j < i 
20 l iL = 1 

where , i = 1 to n 

Next, at step 14, the program solves the 
simultaneous equations for the carrier wave frequency 
under the moment method (wave source of electronic 
25 apparatus and wave source of antenna exist) having the 

mutual impedance obtained by the LDU decomposition at 
step 13 so as to calculate the electric current of the 
carrier wave frequency component flowing through the 
elements of the electronic apparatus segmented at step 4, 
30 solves the simultaneous equations for the upper sideband 

frequency under the moment method (wave source of 
electronic apparatus and wave source of antenna exists) 
having the mutual impedance obtained by the LDU 
decomposition at step 13 so as to calculate the electric 
35 current of the upper sideband frequency component flowing 

through the elements of the electronic apparatus 



- 34 - 



segmented at step 4, and solves the simultaneous 
equations for the lower sideband frequency under the 
moment method having the mutual impedance obtained by the 
LDU decomposition at step 13 so as to calculate the 
5 electric current of the lower sideband frequency 

component flowing through the elements of the electronic 
apparatus segmented at step 4 . 

In solving the simultaneous equations under the 
moment method, the mutual impedance Z is subjected to LDU 
10 decomposition as shown in Fig. 16: 

Z = LDU = LD C L 

From this, the simultaneous equations under the moment 
method are defined, by using the LDU decomposed mutual 
impedance Z , as : 

15 [LD^HI] = [V] 

From this, the simultaneous equations under the moment 
method become equivalent to solving 
[D^Jtl] = [X] and [L][X] = [V] 
The equations can be solved at a high speed since the 

20 matrix is triangularly decomposed. 

From this, by using the LDU decomposed mutual 
impedance common to the carrier wave frequency, the upper 
sideband frequency and the lower sideband frequency for 
solving the three simultaneous equations under the moment 

2 5 method at a high speed (greater number in the case of 
pulse modulation) , the electric current of the carrier 
wave frequency component, the electric current of the 
upper sideband frequency component and the electric 
current of the lower sideband frequency component, 

30 flowing through the elements of the electronic apparatus, 

are calculated at a high speed. 

Next, at step 15, the program uses the electric 
currents found at step 14 in accordance with the above 
method to calculate the voltage generated at a position 

35 between conductor elements specified by the operator, 
then at step 16 judges if the electric current flowing 
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through an element specified by the operator exceeds a 
prescribed threshold and judges if the voltage generated 
at the position between conductor elements calculated at 
step 15 exceeds a prescribed threshold, registers the 
5 results of the above judgement, then returns to step 5 

for the processing of the next carrier wave frequency. 

When the program judges at step 6 that all of the 
frequencies have finished being selected by repeating the 
processings of step 5 to step 16 (except the processing 
10 of step 7), in the above way, the program proceeds to 

step 7 where it outputs the results of the simulation 
obtained at step 10, step 11, step 14 and step 16 to the 

□ output device 3 in accordance with the prescribed manner 

J^ : [ of outputting and ends the processing. 

m 15 In this way, in accordance with the flow of 

sat ; 

s y processing of Fig. 9 and Fig. 10, the simulation program 

}lj 300 simulates the electric current flowing through the 

™' 1 electronic apparatus to which a prescribed intensity of 

p electric field is applied from an antenna and simulates 

tJ 2 0 the voltage generated at the electronic apparatus. 

^ At this time, the simulation program 300 uses the 

\D antenna model registered in the antenna model library 200 

to set the test conditions. The antenna model, as 
explained above, contains the test conditions for 
25 applying the prescribed intensity of an electric field to 

the electronic apparatus, so by using the antenna model, 
the simulation program 300 can immediately apply a 
prescribed intensity of an electric field to the 
electronic apparatus, without trial and error type 
30 processing, and can immediately simulate the effect on 

the electronic apparatus when the prescribed intensity of 
the electric field is applied. 

At step 13 of this flow of processing, the 
simulation program 300 applied LDU decomposition on the 
35 mutual impedance, but it is also possible to apply LU 

decomposition on the mutual impedance. 

That is, as shown in Fig. 17, it is also possible to 



apply LU decomposition to the mutual impedance in 
accordance with the LU decomposition rules for a matrix. 
Here, the following stand for the matrix Dfdij), the 
matrix Lfl^) and the matrix U(u i:j ): 
u ij = z ij- £1 ik u kj 2 

where, E is the sum for k = 1 to (i-1), j = 1 to n, 
i = 1 to j, i^j 

lij = [Zij " Sl ik lkj]/u ji 

where, 2 is the sum for k= 1 to (j-1), i = 1 to n, 
j = 1 to (i-1) , j < i 
lit = 1 

where, i = 1 to n 

The simultaneous equations under the moment method 
are defined, by using the LU decomposed mutual impedance 
Z , as : 

[LU][I] = [V] 

From this, the simultaneous equations under the moment 
method become equivalent to solving 

[U][I] = [X] and [L][X] = [V] 
The equations can be solved at a high speed since the 
matrix is triangularly decomposed. Therefore, the mutual 
impedance may also be subjected to LU decomposition. 

The present invention will be explained in further 
detail below: 

When a radio wave radiated by an antenna is 
modulated by amplitude, if the radio wave is expanded in 
the frequency domain, as shown in Fig. 6A, it is broken 
down into the carrier wave having the frequency f c , the 
upper sideband wave having the frequency (f c +f m ) and the 
lower sideband wave having the frequency (f c -f n ). 

From this, if expressing the electric current 
flowing through an element due to a wave source of the 
frequency (f c -f m ) by [IJ, the electric current flowing 
through an element due to a wave source of the frequency 
f c b Y tne electric current flowing through an 

element due to a wave source of the frequency (f c +f m ) by 
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[I 3 ], the value of a wave source of a frequency (f c -f m ) by 
[Vi], the value of a wave source of a frequency f c by 
[V 2 J, the value of a wave source of a frequency (f c +f m ) by 
[V 3 ], the mutual impedance at a frequency (f c -f m ) by 
5 [Z(f c -f m )], the mutual impedance at a frequency f c by 

[Z(f c )], and the mutual impedance at a frequency (f c +f m ) 
by [Z ( f c +f m ) ] , the electric current flowing through the 
electronic apparatus due to a radio wave radiated by an 
antenna can be found by solving the simultaneous 
10 equations under the moment method of: 

[Zffc-f.JHIx] = [V x ] 

[Z(f c )][I 2 ] = [V 2 ] 

[Z(f c +f m )][I 3 ] = [V 3 ] 

To realize high speed solution of the simultaneous 
15 equations under the moment method, the present invention, 

considering that: 

(f c -f m ) =f c % (f c +f B ) 
adopts the approach of approximation of: 

[Z(f c -f m )] = [Z(f c )3 = [Z(f c +f ffi )] 
2 0 Taking this approach, first, for example, the 

program calculates [Z(f c )]. Of course, it is also 
possible to calculate [Z(f c -f m )], to calculate [ Z ( f c +f n ) ] , 
or to calculate the mutual impedance at another frequency 
close to the carrier wave frequency f c , but it is 
25 preferable to calculate [Z(f c )] since the carrier wave 
frequency f c is positioned at central of the band. 

For example, if calculating [Z(f c )], the electric 
current flowing through the electronic apparatus due to a 
radio wave radiated by an antenna is found by solving the 
30 three simultaneous equations of the moment method of: 

[Z(f c )][Ii] = [V 2 ] 

[Z(fc)][I 2 ] = [V 2 ] 

[Z(f c )J[I 3 ] = [V 3 ] 

In solving the simultaneous equations under the 
35 moment method, if it is possible to ignore the wave 
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source of the electronic apparatus, then only the single 
wave source of the antenna is left as a wave source. Due 
to this, the relation of: 
V 1 :V 2 :V 3 = I^I^Ia 
5 stands, so there is no need to solve the three 

simultaneous equations under the moment method. The 
electric current is found by solving just one among these 
and the electric currents of the remaining frequency 
components are found by proportional operations . 
10 In this case, when the radio wave radiated by an 

antenna is modulated by pulse, there are a plurality of 
upper side frequencies and lower side frequencies, but in 
this case as well just one of the multiple (defined in 
accordance with the above) simultaneous equations under 
15 the moment method is solved and the electric currents of 

the remaining frequency components are found by 
lil proportional operations . 

i=H On the other hand, in solving the above simultaneous 

:~ equations under the moment method, if it is not possible 

y 20 to ignore the wave source of the electronic apparatus and 

therefore the above-mentioned proportional relation does 
=|j not stand between the electric current and the wave 

source, then it is necessary to solve the three 
simultaneous equations under the moment method (in the 
25 case of pulse modulation, a further larger number of 

simultaneous equations). 

In this case as well, since the mutual impedances 
are assumed to be common, it is sufficient to calculate 
the mutual impedance just once, whereby it is possible to 
30 solve the simultaneous equations under the moment method 

at a high speed. 

Further, at this time, LDU decomposition or LU 
decomposition may be applied to the mutual impedance in 
consideration of the fact that the mutual impedance is 
35 common. The processing time required increases due to the 

LDU decomposition or LU decomposition, but use of an LDU 



- 39 - 



decomposed or LU decomposed mutual impedance enables the 
simultaneous equations under the moment method to be 
solved at a high speed. Due to this, when there are two 
or more simultaneous equations under the moment method to 
be solved, the total processing time can be greatly 
reduced. When a radio wave radiated by an antenna is 
modulated by pulse (pulse modulation) , this method is 
extremely effective. 

While not explained in the flow of processing of 
Fig. 9 and Fig. 10, there are cases where the wave source 
of the electronic apparatus can be ignored for some 
simultaneous equations among the plurality of 
simultaneous equations under the moment method. That is, 
in the case of amplitude modulation, there are three 
simultaneous equations under the moment method, while in 
the case of pulse modulation, there are a larger number 
of simultaneous equations under the moment method, but 
sometimes it is possible to ignore the wave source of the 
electronic apparatus for some simultaneous equations 
among them. 

In such a case, while it is necessary to solve the 
simultaneous equations under the moment method for those 
of the wave sources of the electronic apparatus which 
cannot be ignored, it is possible to solve the 
simultaneous equations under the moment method for one of 
those of the wave sources of the electronic apparatus 
which can be ignored and find the electric currents for 
the remaining ones by proportional operations . In this 
case as well, the mutual impedance is calculated only 
once . 

For example, when the basic frequency of a wave 
source of the electronic apparatus is 200 MHz at "f c = 
800 MHz, f m = 1 kHz", the wave source of the electronic 
apparatus for " f c = 800 MHz " cannot be ignored, but the 
wave sources of the electronic apparatus for "f c -f m = 
799.999 MHz " and "f c +f m = 800.001 MHz " can be ignored. 



In this case, the simultaneous equation under the 
moment method for " f c = 800 MHz " is solved and the 
simultaneous equation under the moment method is solved 
for one of "f c -f m - 799.999 MHz " and " f c +f m = 800.001 
MHz", for example "f c -f m = 799.999 MHz " , to calculate the 
electric currents of the frequency component, while the 
electric current component of "f c +f m = 800.001 MHz " is 
found, by a proportional operation, by using the 
calculated electric current of 799.999 MHz, the wave 
source of the antenna at 799.999 MHz and the value of the 
wave source of the antenna at 800.001 MHz. 

In this case, it is necessary to solve two 
simultaneous equations under the moment method in the 
case of amplitude modulation, while it is necessary to 
solve at least two simultaneous equations under the 
moment method in the case of pulse modulation. From this, 
it is preferable to apply LDU decomposition or LU 
decomposition on the mutual impedance as explained above 
and use the decomposed mutual impedance to solve the 
simultaneous equations under the moment method. 

This embodiment assumed the solution of the 
simultaneous equations under the moment method shown in 
Fig. 13 considering only the mutual impedance, but the 
present invention may also be applied as it is to the 
case of solution of the simultaneous equations under the 
moment method shown in Fig. 18 considering the existence 
of a dielectric . 

Solution of the simultaneous equations under the 
moment method shown in Fig. 18 requires calculation of 
not only the mutual impedance Z i;j , but also the mutual 
admittance and the mutual reaction among elements . 
If the simultaneous equations under the moment method 
shown in Fig. 18 are solved, the electric current flowing 
at the surface of the dielectric and the magnetic current 
flowing at the surface of the dielectric are calculated. 

Note that in the equation shown in Fig. 18, I c , n 
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indicates the electric current flowing through metal, I d)I1 
the electric current flowing at the surface of the 
dielectric, Mn the magnetic current flowing at the 
surface of the dielectric, the superscript 0 the value in 
5 air, the superscript d the value in the dielectric, the 

suffix c metal, and the suffix d a dielectric. 

Summarizing the advantageous effects of the 
invention, as explained above, since the apparatus for 
calculating immunity from a radiated electromagnetic 
10 field of one aspect of the present invention divides the 
radio wave radiated by an antenna into a carrier wave, 
upper sideband wave and lower sideband wave and uses the 
□ moment method to simulate the effect of the radio wave 

£\ radiated by an antenna by calculating the mutual 

ID 15 impedance for only one frequency component, using that 

\i mutual impedance to solve the simultaneous equations 

III under the moment method for the carrier wave frequency to 

u * calculate the electric current of the carrier wave 

q frequency component flowing through an electronic 

y~ 20 apparatus due to a radio wave radiated by an antenna, 

solving the simultaneous equations under the moment 
>£? method for the upper sideband frequency to calculate the 

electric current of the upper sideband frequency 
component flowing through the electronic apparatus due to 
25 a radio wave radiated by an antenna, and solving the 

simultaneous equations under the moment method for the 
lower sideband frequency to calculate the electric 
current of the lower sideband frequency component flowing 
through the electronic apparatus due to a radio wave 
30 radiated by an antenna, it is possible to simulate the 

electric current flowing through the electronic apparatus 
at a high speed . 

Further, since the apparatus for calculating 
immunity from a radiated electromagnetic field of another 
35 aspect of the present invention divides a radio wave 
radiated by an antenna into a carrier wave, upper 



sideband wave and lower sideband wave and uses the moment 
method to simulate the effect of a radio wave radiated by 
an antenna by calculating the mutual impedance for only 
one frequency component, using that mutual impedance to 
solve the simultaneous equations under the moment method 
for one of the frequencies, while ignoring the wave 
source of the electronic apparatus, to calculate the 
electric current of that frequency component flowing 
through an electronic apparatus, and calculating the 
electric currents of the remaining frequency components 
by proportional operations, it is possible to simulate 
the electric current flowing through the electronic 
apparatus due to a radio wave radiated by an antenna at a 
high speed. 

Further, since the apparatus for calculating 
immunity from a radiated electromagnetic field of still 
another aspect of the present invention divides a radio 
wave radiated by an antenna into a carrier wave, upper 
sideband wave and lower sideband wave and uses the moment 
method to simulate the effect of a radio wave radiated by 
an antenna by calculating the mutual impedance for only 
one frequency component, using that mutual impedance to 
solve the simultaneous equations under the moment method 
for the frequency overlapping the frequency of a wave 
source of the electronic apparatus to calculate the 
electric current of that frequency component flowing 
through the electronic apparatus, solving the 
simultaneous equations under the moment method for just 
one of the electric currents of the nonoverlapping 
frequency components to calculate the electric current of 
that frequency component flowing through the electronic 
apparatus, and calculating the electric current of the 
remaining frequency components by proportional 
operations, it is possible to simulate the electric 
current flowing through the electronic apparatus due to a 
radio wave radiated by an antenna at a high speed. 

Further, since the apparatus for calculating 
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immunity from a radiated electromagnetic field of still 
another aspect of the present invention uses the moment 
method to simulate the effect of a radio wave radiated by 
an antenna by preparing in advance antenna information 
5 realizing a prescribed intensity of an electric field on 

an electronic apparatus and using that antenna 
information to solve the simultaneous equations under the 
moment method to calculate the electric current flowing 
through the electronic apparatus when there is a 
10 simulation request, it is possible to simulate the 

electric current flowing through the electronic apparatus 
due to the application of the prescribed intensity of an 
i"\ electric field due to a radio wave radiated by an antenna 

at a high speed, 
jri 15 In this way, according to the apparatus for 

n*j calculating immunity from a radiated electromagnetic 

nj field of the present invention, it becomes possible to 

til simulate the electric current flowing through an 

electronic apparatus due to a radio wave radiated by an 
y 20 antenna at a high speed. 

yJ While the invention has been described by reference 

ivi to specific embodiments chosen for purposes of 

illustration, it should be apparent that numerous 
modifications could be made thereto by those skilled in 
2 5 the art without departing from the basic concept and 

scope of the invention. 



